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Abstract 



The minimal supersymmetric standard model (MSSM), is perhaps the leading candi- 
date for new physics beyond the standard model, but it encounters difficulties with 
string gauge unification and in addition does not shed any light on the question of 
fermion masses. We consider a scenario in which the MSSM is valid up to an energy 
scale of ~ 10^^ GeV, but that above this scale the theory is supplemented by extra 
vector-like representations of the gauge group, plus a gauged U{l)x family symmetry. 
In our approach the extra heavy matter above the scale ~ 10^^ GeV is used in two 
different ways: (1) to allow (two-loop) gauge coupling unification at the string scale; 
(2) to mix with quarks, leptons and Higgs fields via spaghetti diagrams and so lead to 
phenomenologically acceptable Yukawa textures. We determine the most economical 
models in which the extra matter can satisfy both constraints simultaneously. We 
then give a general discussion of the infra-red fixed points of such models, pointing 
out the conditions for infra-red stability, then discuss two semi-realistic examples: a 
Higgs mixing model, and a quark mixing model. 



1 Introduction 



The apparent unification of the gauge couphngs at a scale Mqut ~ 10^^ GeV [0 is 
an encouraging feature of the MSSM, whose direct experimental verification so far 
remains out of reach. However the MSSM does not address the question of fermion 
masses, and in any case it is unlikely that the MSSM can survive unchanged above 
Mgut ~ 10^^ GeV since the gauge couplings which converge on this scale begin to 
diverge above it, and are quite unequal at the string scale Mx ~ 5 x 10^'' GeV, 
even taking into account higher Kac-Moody levels and string threshold effects [Q. 
The traditional approach is to embed the MSSM in some supersymmetric grand 
unified theory (SUSY GUT) but such an approach presents many theoretical and 
phenomenological challenges [0, ^, and we shall not consider it further here. 

String models predict that the three gauge couplings of the MSSM are directly 
related to each other at the string scale Mx [§!• Weakly coupled string theories give 
the relation |^ 

Mx = 5.27 X lO^V GeV, (1) 

where gx is the unified gauge coupling at the string scale Mx- If the MSSM (and 
nothing else) persists right up to the string scale Mx such theories are in conflict 
with low energy measurements of the gauge couplings. From this data, one can show 
that gauge couplings cross at ~ 10^^ GeV, and significantly diverge at the string 
scale ~ 5 X 10^'' GeV. However the situation is in fact not so clear cut since the 
U{1)y hypercharge gauge coupling has an undetermined normalisation factor ki > 1 
(where for example ki = 5/3 is the usual GUT normalisation) which may be set 
to be a phenomenologically desired value by the choice of a particular string 
model. However the simplest string theories (e.g. heterotic string with any standard 
compactification) predict equal gauge couplings for the other two observable sector 
gauge groups g2 = at the string scale Mx, which would require a rather large 
correction in order to account for ^^(m^) |^. In fact, a recent analysis p|, |10[ 
concludes that string threshold effects are insufficient by themselves to resolve the 
experimental discrepancy. The analysis also concludes that light SUSY thresholds 
and two-loop corrections cannot resolve the problem, even when acting together. In 
order to allow the gauge couplings to unify at the string scale it has been suggested 



111 that additional heavy exotic matter in vector-like representations should be added 
to the MSSM at some intermediate scale or scales Mj < Mx, leading to the so called 
MSSM-I-X models. A detailed unification analysis of such models was performed by 
Martin and Ramond (MR) for example. In a previous paper |TB| we performed 
a general unification analysis of MSSM-I-X models, focusing on the infra-red fixed 
point properties of the top quark mass prediction, using similar techniques to those 
proposed for the MSSM and GUTs [|1^, |16[. The main result was that the top 
quark mass tends to be heavier than in the MSSM, and closer to its quasi-fixed point 
in these models. 

This leads us to the question of fermion masses in the MSSM. The fermion 
mass problem in the MSSM arises from the presence of many unknown dimensionless 
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Yukawa couplings. A possible solution to this long-standing problem is that the 
MSSM could be embedded in some more predictive high energy theory with fewer 
Yukawa couplings. Alternatively one may attempt to relate the Yukawa couplings to 
the gauge coupling, as in the Pendleton-Ross infra-red stable fixed point (IRSFP) for 
the top quark Yukawa coupling ||T3[, or the quasi-fixed point of Hill To obtain 



IRSFPs it is necessary that the unknown dimensionless Yukawa couplings are of the 
same order as the gauge coupling at high energy. Thus at first sight this approach 
would seem to be inapplicable to the small Yukawa couplings of the standard model. 
However string theory yields large Yukawa couplings of order the gauge coupling, and 
this could suggest that the small dimensionless Yukawa couplings be reinterpreted as 
large Yukawa couplings multiplied by a ratio of mass scales. In such a scenario fixed 
points may once again be applicable, as pointed out recently by Ross [|18]. In this 
case the high energy scale would be the string scale and the infra-red region would 
be near the GUT scale. Thus IRSFPs may provide an unexpected resolution of the 
fermion mass problem. 

In the present paper we consider a bottom-up string-inspired approach in which 
the MSSM is valid up to an energy scale of M/ ~ 10^^ GeV. Above this scale the 
theory is supplemented by extra vector-like representations of the gauge group, plus a 
gauged U{l)x family symmetry, which turns out to be pseudo-anomalous as discussed 
in the next section. Although we do not derive our results from a string model, both 
the presence of extra U{1) gauge groups (one of which is pseudo-anomalous and 
may be identified with the gauged family symmetry U{l)x), and extra vector-like 
matter, are typical of recent string compactificationsQ. The basic idea of our present 
approach is that the extra heavy matter above the scale ~ 10^^ GeV may be used in 
two different ways: (1) to allow (two-loop) gauge coupling unification at the string 
scale; (2) to mix with quarks, leptons and Higgs fields via spaghetti diagrams and 
so lead to phenomenologically acceptable Yukawa textures. We emphasise that in 
our approach the operators required for Yukawa textures are generated from the 
dynamics of the effective field theory beneath the string scale rather than from the 
string theory itself. The overall philosophy of the approach thus far is to explain 
the data on fermion masses and mixing angles qualitatively. Unfortunately, once 
the extra fields have been added there are many extra free parameters than data 
points. In order to increase the predictivity of the class of models we give a detailed 
discussion of the IRSFP properties of the models. As indicated above, the IRSFP 
approach has the potential to predict all of the masses and mixings in terms of about 
two free parameters, although as we shall see there are difficulties in practice with 
implementing such a scheme. 



In this paper we essentially build on the work of ref.|]T8[, where a specific model 
with a gauge U{l)x family symmetry was introduced, and the extra vector-like matter 
that was introduced consisted of extra Higgs doublets which mixed with the ordinary 



^The relation between the fixed point and quasi- fixed point has been fuUy elucidated by Lan- 
zagorta and Ross . The conditions for the stability of fixed points for any number of dimensionless 



Yukawa couplings in general supersymmetric models has also recently been examined |17 
^For more discussion about this point see ref.p^ 
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Higgs doublets. Our present work is more general and differs in several respects. For 



example the question of gauge coupling unification was addressed in ref . [|18[ by adding 
complete SU{5) 5 © 5 representations (which contain the additional Higgs states 
required for the mixings) to the MSSM theory with masses just below the unification 
scale. These have no relative effect on the running of the three gauge couplings to 
one loop order, however at two-loop order it was claimed that the unification scale is 
raised. In our two loop analysis we explicitly find that such a mechanism does not 
increase the unification scale sufficiently, so we abandon this possibility here. Instead 
we perform a general two-loop analysis of string gauge unification, including heavy 
matter at a scale M/ which consists of MSSM vector representations which do not 
form complete SU{5) representations, i.e. we allow split SU (5) representations. In this 



respect our approach resembles that in refs. [|ri|, |T2|, |T3[ , however we emphasise that in 
our present approach we impose the additional constraint that the intermediate mass 
scale Mf must be very close to the string scale Mx because, as we shall see in the 
next section, the U{l)x gauge symmetry is pseudo-anomalous and must be broken 
not far below the scale Mx- In practice this implies that Mj ~ 10^^ GeV, as stated 
earlier. Another difference between our work here and ref.[T^ is that we allow more 
general mixing possibilities along the Higgs, quark and lepton lines of the spaghetti 



diagrams, rather than just the Higgs mixing considered in ref.[jT8[. Finally we discuss 
the IRSFPs in this more general framework, rather than just the particular Higgs 
mixing model in ref. [ 



We emphasise that in our approach the extra states required for unification are 
also used for spaghetti mixing. This economical double use of the extra heavy matter 
is the main new idea of the present paper, and one of our main results is to tabu- 
late the minimal amount of extra vector-like matter which allows both string gauge 
unification (subject to the high Mj constraint) and phenomenologically acceptable 
Yukawa textures (see Tables ^ and §.) 

The layout of the remainder of this paper is as follows. In section 2 we show how 
Yukawa textures may be generated from a broken U{l)x gauged family symmetry 
where the desired operators are generated from the effective field theory beneath the 
string scale via spaghetti diagrams. In section 3 we discuss a two-loop analysis of 
string gauge unification, where the extra states are consistent with the requirements 
of the previous section. In section 4 we give a general discussion of the infra-red 
fixed point nature of such models, then discuss two semi-realistic examples. Section 5 
concludes the paper. Appendix 1 summarises the running of the gauge couplings to 
two loops in models with matter additional to the MSSM, appendix 2 lists the wave- 
function renormalisations for the general model discussed in section 4, appendix 3 
details the infra-red fixed point of the Higgs mixing model, and appendix 4 details 
the infra-red fixed point of a quark-line mixing model. 
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2 Yukawa Textures from U{l)x Family Symmetry, 
Vector Representations and Spaghetti Diagrams 



The pattern of quark and lepton masses and quark mixing angles has for a long time 
been a subject of fascination for particle physicists. In terms of the standard model, 
this pattern arises from three by three complex Yukawa matrices (54 real parame- 
ters) which result in nine real eigenvalues plus four real mixing parameters (13 real 
quantities) which can be measured experimentally. In recent years the quark and 
lepton masses and mixing angles have been measured with increasing precision, and 
this trend is likely to continue in the future as lattice QCD calculations provide in- 
creasingly accurate estimates and B-factories come on-line. Theoretical progress is 
less certain, although there has been a steady input of theoretical ideas over the years 
and in recent times there is an explosion of activity in the area of supersymmetric 
unified models. This approach presumes that at very high energies close to the unifi- 
cation scale, the Yukawa matrices exhibit a degree of simplicity, with simple relations 
at high energy corrected by the effects of renormalisation group (RG) running down 
to low energy. For example, the idea of bottom-tau, top-bottom-tau and even top- 
bottom-tau-tau-neutrino Yukawa unification have received a good deal of attention 
recently (see for example ref . ||20|| . ) 



These successes with the third family relations are not immediately generalisable 
to the lighter families. For the remainder of the Yukawa matrices, additional ideas 
are required in order to understand the rest of the spectrum. One such idea is that 
of texture zeroes: the idea that the Yukawa matrices at the unification scale are 



rather sparse; for example the Fritzsch ansatz ||21[|. Although the Fritzsch texture 



does not work for supersymmetric unified models, there are other textures which do. 



for example the Georgi-Jarlskog (GJ) texture p2| for the down-type quark and lepton 
matrices: 

/ Ai2 \ / Ai2 \ 

A^ = A21 -3A22 , A^ = A21 A22 . (2) 
V A33 ; V ^33 / 

After diagonalisation this leads to A^ = A;,, A^ = SXg, Ag = Ad/3 at the scale Mqut 
which result in (approximately) successful predictions at low energy. Actually the 
factor of 3 in the 22 element above arises from group theory: it is a Clebsch factor 
coming from the choice of Higgs fields coupling to this element. 

It is observed that if we choose the upper two by two block of the GJ texture 
to be symmetric, A12 = A21, and if we can disregard contributions from the up-type 
quark matrix, then we also have the successful mixing angle prediction 



Vus = \f\J\s. (3) 



The data therefore supports the idea of symmetric matrices, and a texture zero in the 
11 position. Motivated by the desire for maximal predictivity, Ramond, Roberts and 
Ross (RRR) p3| have made a survey of possible symmetric textures which are both 
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consistent with data and involve the maximum number of texture zeroes. Assuming 
GJ relations for the leptons, RRR tabulated five possible solutions for the up-type 
and down-type Yukawa matrices. 

Having identified successful textures^, the obvious questions are: what is the 
origin of the texture zeroes? and: what is the origin of the hierarchies (powers of 
the expansion parameter e)? A natural answer to both these questions was provided 
early on by Froggatt and Nielsen (FN) |2^. A specific realisation of the FN idea was 
provided by Ibanez and Ross (IR) based on the MSSM extended by a gauged 
family U{l)x symmetry with 6 and 6 singlet fields with opposite X charges, plus 
new heavy Higgs fields in vector representations^. Anomaly cancellation occurs via a 
Green-Schwarz-Witten (GSW) mechanism, and the U{l)x symmetry is broken not far 
below the string scale The idea is that the U{l)x family symmetry only allows 
the third family to receive a renormalisable Yukawa coupling but when the family 
symmetry is broken at a scale not far below the string scale other families receive 
suppressed effective Yukawa couplings. The suppression factors are proportional to 
powers of the vacuum expectation values (VEVs) of 6 fields which are MSSM singlets 
but carry U{l)x charges and are responsible for breaking the family symmetry. The 
non-renormalisable operators that give the small effective Yukawa couplings are scaled 
by heavier mass scales Mj identified as the masses of new heavy vector representations 
which also carry U{l)x charges. IR envisaged a series of heavy Higgs doublets of 
mass Mj with differing U{l)x charges which couple to the lighter families via sizable 
Yukawa couplings that respect the family symmetry. The heavy Higgs doublets also 
couple to the MSSM Higgs doublets via 6 fields and this results in suppressed effective 
Yukawa couplings when the family symmetry is broken. 

By making certain symmetric charge assignments (see later), IR showed that the 
RRR texture solution 2 could be approximately reproduced. To be specific, for a 
certain choice of U{l)x charge assignments, IR generated Yukawa matrices of the 
form: 



(4) 



These are symmetric in the expansion parameters e and e, which are regarded as inde- 
pendent parameters. This provides a neat and predictive framework, however there 
are some open issues. Although the order of the entries is fixed by the expansion 
parameters, there are additional parameters of order unity multiplying each entry, 
making precise predictions difficult. A way to address the problem of the unknown 
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coefficients has been proposed in |]I8[ where it has been shown that the various coeffi- 
cients may arise as a result of the infra-red fixed-point structure of the theory beyond 
the Standard Model. The idea behind this approach is that since there are no small 
Yukawa couplings one may hope to determine all the Yukawa couplings by the use 



■'Over the recent years, there has been an extensive study of fcrmion mass matrices with zero 
textures Q. 

^The generahsation to include neutrino masses is straightforward [Em. 
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of infra-red fixed points along similar lines to the top quark Yukawa coupling deter- 
mination. This attractive idea thus allows the determination of otherwise unknown 
parameters in the model, simply by the dynamics of the renormalisation group (RG) 
flow of the model, as discussed in Section 4. 

Having given an overview of the approach, we now discuss in a little more detail 
how the U{l)x gauge symmetry is introduced. The way this is achieved is well 
documented and here we only sketch the main results. The quark and lepton 
multiplets are assigned family dependent (FD) charges as shown in Table |I]. 





Q, Uf Df U Ef Hi H2 


U{1)fd 


ai ai ai ai ai -la^ -la^ 



Table 1: U{\)pi) charges assuming symmetric textures. 

The need to preserve SU{2)l invariance requires left-handed up and down quarks 
(leptons) to have the same charge. This, plus the additional requirement of symmet- 
ric matrices, indicates that all quarks (leptons) of the same i-th generation transform 
with the same charge a^. It is further assumed that quarks and leptons of the same 
family have the same charge (this additional assumption was made by Ross ||18|| ). 



The full anomaly free Abelian group involves an additional family independent com- 
ponent, U{1)fi, and with this freedom U{1)fd is made traceless without any loss of 
generality^. Thus we set ai = —(02 + ^3). 

Making the above assumptions all charge and mass matrices have the same struc- 
ture under the U{1)f£) symmetry. The FD charge matrix involving two quark or 
lepton fields and a Higgs is of the form 

/ — 2^2 — 4^3 — 3a;3 —0:2 — 2^3 \ 

— 3^3 2(^2 — «3) ^2 — «3 (5) 
\ —a2 — 2^3 a2 — as / 

Acceptable textures are obtained for 

^^=1 (6) 

a2 — as 

or 

"2 = 2^3 (7) 

However these are only the family dependent charges. The total U{l)x charges 
are given by 

f/(l)x = f/(l)FD + f/(l)F/ (8) 

where 

Uil)FI = Uil)T + Uil)F (9) 

where the corresponding FI charges are denoted by t,f and the resulting U{l)x 
charges are given in Table |^. 

^Since we assume that the 33 operator is renormalisable, the relaxation of the tracelessness 
condition does not change the charge matrix since any additional FI charges can always be absorbed 
into the Higgs charges Hi 2- 
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u 




H2 


Uil)x 


ai + t ai + t 


ai + f 


ai + f 


ai + t -2^3 - (/ + t) 


-2a-i - 2t 



Table 2: U{l)x charges assuming symmetric textures. 



At this point we should be clear about the question of anomalies for the various 
?7(l)'s which were introduced into this model [^. We began by introducing a.U{l)FD 



symmetry, which was arranged to be traceless without loss of generality. Although 
the tracelessness is sufficient to guarantee that the quarks and leptons contribute a 
zero anomaly, the phenomenologically required Yukawa texture implies that the Higgs 
doublets must carry a U{1)fd charge (the above conditions cannot be satisfied by 
as = 0) and consequently the Higgs contribute a non-zero SU{2YU{l)p£) anomaly. 
In order to improve matters, we returned to the U{1)fi symmetry which is a family 
independent symmetry, identified as the sum of U{1)f and U{1)t- Although there is 
no explicit SU{5) unification, the components of the "5" rep are assigned a charge / 
under U{1)f, while the components of the "10" are assigned a charge t under U{1)t, 
and the Higgs doublets are assigned the charges which guarantee that the Yukawa 
charge matrix remains unchanged under the total U{l)x symmetry defined above as 
the sum of the U{1)fd and U{1)fi symmetries. Although, following ref.|^, we have 



introduced several ?7(l)'s, it should be remembered that at the end of the day only 
one of them (the total U{l)x) is regarded as the physical gauged family symmetry, 
and the others are for construction purposes only. The essential point to note is that 
the choice of FI charges in Table ^ is the most general choice consistent with anomaly 
cancellation amongst the quarks and leptons via the Green-Schwarz mechanism |^ 
which is based on assigning SU{5) multiplets equal charges. The physics of the GS 
mechanism in 4D is that one cancels the anomaly of the f/(l)x by an appropriate shift 
of the axion present in the dilaton multiplet of 4D strings. For the GS mechanism to 
be possible, the coefficients of the mixed anomalies of the U{l)x with SU{3), SU{2) 
and U{1)y have to be in the ratio of the Kac-Moody levels : k2 : ki, and hence 
the GS mechanism works if the anomalies are in the ratio 1:1:5/3, and this is in 
fact guaranteed by the choice of quark and lepton charges in Table ^ This of course 
still leaves the question of the contribution to the anomalies from the Higgs doublets. 
In order to allow anomaly cancellation amongst the Higgs doublets, and permit a 
renormalisable /i mass term between the Higgs doublets we must arrange that the 
Higgs carry zero X charge, which in turn implies that the third family has zero X 
charge, and this is accomplished by choosing: 

f = t = -as, (10) 

Putting 03 = 1, where the U{l)x symmetry is broken by the VEVs of MSSM singlet 
fields 6 and 6 with U{l)x charges -1 and +1 respectively implies the following simple 
X charges for the three families: 

1"* family: X = -4 

2"'^ family: X = +1 . (11) 
3'''^ family: X = 
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Since the Higgs charges are zero the charge matrix is simply the matrix given by the 
sum of the quark (or lepton) charges from each family: 




(12) 



In any mass term, the sum of the charge of the fields in that term must be zero to 
preserve U(l)x- Thus from Eq.|l^, the operators that could possibly generate the U 
quark mass matrix are: 

Q2U^H2{9f, Q3U^H2{9)\ QiU^H2{9)\ 
QiU1H2{ef, Q2U^H2{e)\ QiU^H2{ef. 



(13) 



Because each 6 or 6 field corresponds to a suppression factor of e (where e =< 6 > 
/Ml =< 6 > /Mj ) the corresponding texture is of order 



/ 




(14) 



We now turn to the question of the origin of the non-renormalisable operators in 
Eg .[Tsl. A natural answer to this question was provided by Froggatt and Nielsen (FN) 
|25| . The basic idea involves some new heavy matter of mass^ M/ which are in vector 



representations of the MSSM gauge group and which carry charges under the family 
group U{l)x- The vector-like matter couples to ordinary matter (quarks, leptons, 
Higgs) via the MSSM singlet fields 9 and 9 leading to "spaghetti-like" tree-level di- 
agrams. The spaghetti diagrams yield the effective non-renormalisable operators in 



Eq.l3 



An explicit realisation of this mechanism was discussed by IR^3 subsequently 
by Ross [|18] , based on the heavy vector-like matter corresponding to additional Higgs 
doublets which could mix with the MSSM doublets Hi, H2 via the spaghetti diagrams. 



Thus the following Higgs were introduced |18 



o-(-l) frW 7t(-2) fr(2) rr{3) o-(-3) ^.(4) f.(-4) ^-(S) fr{-S 
-"1,2 5-"l,2)-"l,2 ) -"1,2 5 -"1,2 ) -"1,2 )-"l,2?-"l,2 5-"l,2)-"l,2 

(x) 



(15) 



where the U{l)x charges are given in parentheses, and Hi 2 have hypercharges Y/2 = 
— 1/2,1/2, and h[~2^^ have hypercharges Y/2 = 1/2,-1/2, respectively. The idea is 

(x) 

that the Higgs Hi 2 have direct couplings to the lighter families and mix with the 



^We will identify the new heavy matter required to generate the non-renormalisable operators 
with the new heavy matter required to ensure string gauge unification so we take the mass scale to 
be Mj. 
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H 



(0) 



\/ 



Figure 1: Renormalisable 33 operator. 



MSSM Higgs Hi 2 via singlet 9 fields. Thus the renormalisable Higgs terms (where 
we have neglected the Yukawa couplings) are : 



( Qi Q2 Qz 



rriS) 
-"1,2 
^(3) 
-"1,2 

rrW 

-"1,2 



rj{3) 
-"1,2 

rr{-2) 

-"1,2 

o-(-l) 

-"1,2 



Hi'} \ 

Hl,2 j 



1) 



(16) 



For example, the renormalisable 33 operator is depicted in Fig.|l|. 

Such direct Higgs couplings, allowed by the U{\)x symmetry, combined with Higgs 
mixing via singlet field insertions, lead to the effective non-renormalisable operators 
in Eq.|l2|. For example the spaghetti diagram responsible for the 32 quark mixing 
term is illustrated in Fig.^. It is clear that such diagrams generate all the operators 
in Eq.|l^. By drawing such diagrams it becomes clear that additional heavy vector 
Higgs are required beyond those in the direct coupling matrix in Eq.p!6|. It is easy to 
see that all Higgs charges in integer steps between 8 and -2 are required if all elements 
of the mixing matrix are to be non-zero. For example Cabibbo mixing requires the 
additional Higgs with charges 2 and 1 (plus their conjugates) so that the Higgs of 
charge 3 can step down to the Higgs of zero charge via three 9 field insertions as 
shown in Fig.^. 

The full list of Higgs required for this scenario is larger than assumed in Eq.|l5| 
and is displayed below: 



rr{8) o-(-8) rr(7) fT{~^) tt{^) fr(-6) rr(5) fr(-5) rr(4) fr(-4) 
-"1,2!-"1,2 ?-"l,2i-"l,2 )-"l,2)-"l,2 )-"l,2)-"l,2 i-"l,2?-"l,2 ! 

H^^^ H 
-"1,2 5 -" 



(-3) tt{2) fr{-2) rr(l) fr(-l) rr(-l) fr(l) tt{-2) fr(2) 
1,2 5 -"1,25 -"1,2 5 -"1,25 -"1,2 5 -"1,2 5 -"1,25 -"1,2 5 -"1,2 



(17) 



Since complete Higgs mixing requires 20 Higgs vector representations, rather than 
10 as assumed on the basis of the direct Higgs couplings, it is natural to try and 
look for a more economical alternative. This provides a motivation to study quark 
and lepton mixing in addition to Higgs mixing, as a means of generating the desired 
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Figure 2: Spaghetti diagram for 32 mixing 



non-renormalisable operators. For example in Fig.^ we show a spaghetti diagram 
which can generate 32 mixing along the line. Notice that we have added an 
intermediate quark f/*-"^ with the same quantum numbers under the gauge symmetry 
of the model as U^. In general, both will mix with U^^^ through a heavy mass term 
and we may always rotate the definition of the fields such that one of the linear 
combinations is massless. We identify this combination with the (conjugate) right 
handed top superfield of the MSSM and the other with the heavy field involved in 
the mass suppression of the spaghetti diagrams]]. This means that vector quark and 
lepton states must be added, one chiral partner of which has the same X charge as a 
MSSM state. In Fig.^ we show how 32 mixing can be generated by mixing along the 
Q line. 

Similar diagrams can be drawn for the other mixings, and by consideration of such 
diagrams we can see that vector fields with X charges from -4 to 4 in integer steps 
are required. For example complete mixing along the Q line requires: 

g«, g(-^\ Q(^\ Q^-'\ Q^'\Q^-'\Q^'\Q^-'\ (18) 

corresponding to 9 vector Q + Q representations. Also, mixing along the line 
requires 9 vector + representations; mixing along the -D^ line requires 9 vector 
D'^ + D'^ representations; mixing along the L line requires 9 vector L + L represen- 
tations; and mixing along the E'^ line requires 9 vector E'^ + E^ representations. In 

''Wherever two fields with identical quantum numbers are present, one of which is a state of 
the MSSM, we will assume that this mixing has already been accounted for. The states labeled as 
MSSM states will theforc be defined as the massless linear combination. 
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ure 3: Spaghetti diagram for Cabibbo (12) mixing 
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Figure 4: Spaghetti diagram for 32 mixing along the U hne 
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Figure 5: Spaghetti diagram for 32 mixing along the Q hne 



each case the X charges run from -4 to +4 in integer steps (plus the opposite charges 
for the conjugate states), as in EqJT^. 

Now that we have allowed quark and lepton mixing (along the doublet and/or 
singlet lines) as well as Higgs mixing, many possibilities present themselves, corre- 
sponding to different combinations of each type of mixing. In the next section we 
shall use the constraint of string gauge unification to help to discriminate between 
the different possibilities. Here we shall make some general observations about the 
model building. 

Let us begin with the U mass matrix. We can envisage a scenario in which we 
have the ^ Q + Q fields listed in Eq.|l^, plus 9 vector U'' + U'' representations. In 
addition to these we also have the three chiral families with the X charges -4,1,0 as 
discussed above. To be more general we must also consider additional Higgs vector 
representations. We have seen that a possible Higgs mixing scenario requires 10 
H2 + H2 plus 10 Hi + Hi extra Higgs with H charges from 8 to -2. However we 
may also wish to consider Higgs which couple any vector quark field to any other 
vector quark field, in which case the Higgs charges must range from 8 to -8. If for 
the moment we ignore the MSSM fields, but include all of the extra vector fields 
mentioned above then we have a 9 x 9 matrix of Higgs couplings to quark fields: 
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(19) 



There is of course a similar matrix of couplings involving the conjugate fields. Let us 
now introduce the three families of MSSM chiral fields shown below into the matrix 



in Eq.l9 



-4) rr(l) tAO) 



(20) 



where the family index is indicated by a subscript and the X charge is indicated by a 
superscript. The matrix now becomes a 12 x 12 matrix, and the Higgs content stays 
the same. We can now consider all possible ways in which mass mixing between the 
MSSM quarks can occur. In general for ij mixing we require a spaghetti diagram 
with the external lines consisting of Qi, Uj, ■ In addition we are allowed to hang 
any amount of 6 and 6 spaghetti along any of the three lines Q,U,H2 in order to 
achieve the mixing where the minimum amount of spaghetti corresponds to the lead- 
ing non-renormalisable operator. To take a trivial example the 33 operator Q^U^ut 
is achieved directly at tree level without any spaghetti. At the other extreme the 11 
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operator Qi ^^f/} is clearly forbidden at tree level by the X symmetry, with 

the allowed non-renormalisable operator being (5i~^^f/i~^^-ff2°''(^''^'')^- The required 8 
pieces of 6 spaghetti can be hung along any of the three hues Q,U,H2 depending on 
the vector fields and charges which are assumed. For example in the Higgs mixing 
scenario of IR there is a tree level Higgs coupling Q[~^^ Ui~^^ H2^^ and the MSSM Higgs 
iJ2°^ is achieved by stepping down the Higgs charge along the first row of the matrix 



in Eq.l9 with the Higgs charge decreasing by one unit after each 9 field insertion: 



^(8) ei^ ^(7) ^(6) ^(5) ^(4) ei^ ^(3) ei^ ^(2) ei^ ^(i) 0^ ^(o) ^21) 

with all the spaghetti mixing along the Higgs line. With the additional vector Q and 
fields considered above there are alternative ways in which the spaghetti mixing 
can take place. For example we could begin from the non-renormalisable operator 
jj(4)^^o) g^^p down the the U line to reach the desired f/^"^'' field: 

fj(i) (jW (jW (jW fjW (jW fj(i) (jW , 

f/(4) % JjiS) % jj{2) \ \ jj{0) % jji^l) \ jj{-2) \ jj{-3) \ ^{-i) ^22) 

The resulting spaghetti diagram now has all the mixing along the U line, but is of 
the same order as the previous diagram. We could repeat this starting instead from 
the non-renormalisable operator Q'^'^^u[ ^''H'f''^ and step down the the Q line to reach 
the desired Qi"^'' field: 

g(4) % g(3) % g(2) % g(l) % g(0) % g(-l) % g(-2) % g(-3) % g(-4) (23) 

The resulting spaghetti diagram now has all the mixing along the Q line, but is of 
the same order as the previous diagram. There are of course many other possibilities 
which involve a combination of the three types of mixing, and all these possibilities 
will lead to non-renormalisable operators of the same order. For example suppose we 
again wish to generate the 11 operator Q'\~^^u[~^'^ H'2\o'^'^''Y starting from the tree 
level operator Q^~'^'>U''~'^^ Now in order to achieve this we must have all three 
types of mixing simultaneously: 

^(4) ^(3) ^(2) ^(1) ^(0) (24) 

(25) 
(26) 

Again the 11 operator is eighth order, but now there are four pieces of spaghetti 
from the Higgs line, two from the U line and two from the Q line. There are clearly 
many other possible ways of achieving 11 ?7 mixing. The discussion of the other U 
mixings is similar. Finally the discussion of the D and E mixing matrices follows a 
similar pattern. The most general model clearly involves 9 vector copies of each of 
{Q + Q),iU + U),{D + D),{L + L),{E + E), where we denote the number of vector 
copies as nQ,nu,nD,nL,nE respectively, plus 16 vector copies of each of the Higgs 
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fields {H2 + H2), {Hi + Hi), where we denote the number of vector copies as nHi,nH2, 
respectively. 

This is clearly not the most economical model. For example the Ross[0 model 
is based on no extra vector quarks and leptons and 10 vector copies of each of {H2 + 
H2), [Hi + Hi). In the next section we shall impose the constraint of string gauge 
unification in order to try to determine a more economical model. However it is clear 
from the discussion of this section, that if too few extra vector states are allowed, then 
the required mass mixing will not be achievable. Therefore we seek the minimal model 
which is consistent with the constraints of spaghetti mixing discussed here, and string 
gauge unification. Since mixing can be achieved by a combination of mixing along 
the three different lines of the spaghetti diagram, in the next section we shall impose 
the following conservative lower limits on the minimal numbers of vector copies of 
fields. From U mixing we require: 

nQ + n2 + nu> 8. (27) 

From D mixing we require: 

nQ + n2 + nD> 8. (28) 

From E mixing we require: 

n2 + nE>8. (29) 
where for convenience we have defined the total number of doublets as 

712 = riH^ + + riL- (30) 

In addition we shall allow for exotic superfields known as "sextons" S which are colour 
triplets, electroweak singlets and have hypercharge 1/6. The sextons occur together 
with their vector conjugate {S + S) and we denote the number of vector copies of 
sextons as ns- These superfield representations are present in the massless spectrum 
of some string models [ ^ . 



3 String gauge unification analysis 

We now describe the numerical constraints placed upon the models to ensure that 
they provide agreement with the low energy data, string scale gauge unification and 
are compatible with models of fermion mass and mixing. Since we are hoping to 
eventually embed the model into a free-fermionic Kac-Moody level 1 string model, 
we must make sure that the gauge couplings obey the constraint of gauge unification 
at the string scale as in Eq.|I|. Another constraint comes from the agreement with the 
empirically determined values of the gauge couplings at low energy scales |^ . The 
data used is 

a{MzY^ = 127.90 ±0.09 

sin^^^ = 0.2315 ±0.0002 
as{Mz) = 0.118 ±0.003, (31) 
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where the numbers quoted are those derived in the MS renormahsation scheme from 
experiments. In what follows, we shall assume central values for a{Mz)~^ Ow 
because their errors are comparatively small. The third constraint comes from the 
fact that we are expecting to use the intermediate matter as the heavy fields in 
a family U(l)x model of fermion masses. As previously demonstrated these 



models require Mi to be within a couple of orders of magnitude of Mx to make 
the GS gauge anomaly cancellation mechanism work consistently. After some other 
filtering of models, as described below, the condition imposed on any successful model 
will be 

Mi/Mx > 1/40. (32) 

Previous models of U(l)x family symmetry require a GUT type normalisation of Y 
from the anomaly cancellation conditions. Therefore, our condition upon the unifica- 
tion of gauge couplings will be 

g,{Mx) ^ yf = 92{Mx) = gsiMx), (33) 

corresponding to a Kac- Moody level 1 string model with ki = 5/3. Finally, bearing 
in mind the long-term view of requiring the intermediate sector to mediate masses 
and mixings to all of the SM fermions, we impose Eqs.^^^-|29| We then search through 



all of the models satisfying Eqs. |27| - pp| for kq, ?t,2, riu, nj:,^ ue < 10 in order to find 
the models with less field content. 

We now describe the systematic procedure to determine the models that pass 
the constraints^ given by Eqs. p7| - p^ , pl| - |33| . Because of computer time constraints, we 
were not able to determine the predictions to two-loop order of every model satisfying 
Eqs. p7| - P^ . The following procedure was therefore adopted: the predictions for every 
different choice of intermediate field content were determined to one loop order as 



in [T^l for ki = 5/3. If the one-loop predictions did not satisfy certain constraints to 
be described shortly, the models were discarded. For any models passing the previous 
"cut", the predictions were obtained at two-loop order. 

The first one-loop filtering procedure was as follows: once a model has been se- 
lected by a a particular choice of intermediate matter, a guess (Mx') of the string 
unification scale was made. The condition ai{Mx') = a2{Mx') yields a value of Mj 
consistent with unification at Mx by solving the one-loop RGEs for ai and a2 in the 
MS scheme to obtain 

InMj = (2TT{a^\Mz)-a^\Mz))-—\nMz + —\nmt + -\nMsusY 

56 fi \ 3 

-(3nQ + ^2 - - - -Yt) InMx'j /{Yt- - Suq - ris), (34) 

where Yt = J^iiXif^Y, i runs over all of the intermediate states and Yi denotes 
the hypercharge of the intermediate state i. Throughout the numerical analysis, the 
entire SUSY spectrum was assumed to be at an effective scale equal to Msusy- Any 

®In fact, we will allow the prediction of as{Mz) to be within 2a of the value quoted in Eq.pl|. 
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threshold effects were taken into account by a step function approximation. With M/ 
and Mx values consistent with gauge unification at a scale Mx', we could calculate 
Mx consistent with string scale unification in Eq.0 by finding 

^2 [Mx) = [Mz) + — - In — + — - In 7^ In ■ 



127C Mz 127r mt 27r M, 



SUSY 



(35) 



2n Mi 

and substituting it into Eq.|l]. To obtain values of Mx,Mj consistent both with 
ai{Mx) = a2{Mx) and Eq.[|, we now substitute Mx' with Mx and iterate the above 
procedure until Mx = Mx' is satisfied to some required accuracy. This yields a 
prediction for a^i^Mz) by using a^i^Mx) = O-g where ac is the string scale unified 
gauge structure constant, is then run to low energies using the one loop RGEs, 

23 1 2 

a^\Mz) = aG{Mx) + — hiMz-—\nmt--\nMsusY- 

OTT OTT TT 

2n^ + + ^^^^ _ 3 - 2ng - - ^^^^ ^^^^ 
2tx 2tt 

We have set Msusy = f^^t-, whichQ we take to be 166 GeV, corresponding to a top 
quark pole mass of 180 GeV for central values of as{Mz)-, as in ref. We will 

return later to the effect of the empirical errors upon the inputs. We now require each 
model to pass the cuts as{Mz) < 0.124 and Mi/Mx > 1/100 to be worthy of the 
two-loop analysis. Note that these constraints are purposefully less severe than the 



ones in Eqs. |31| , |32| because we do not want to discard models in which the imprecise 
one-loop predictions do not pass the more restrictive constraints, but in which the 
two-loop predictions pass. 

Having attained a list of all models that passed the initial cuts, the two-loop pre- 
dictions were then attained. At the two-loop level, the third family Yukawa couplings 
all effect the running of the gauge couplings and therefore the predictions of gauge 
unification. To a good approximation, the other Yukawa couplings of the MSSM 
have a negligible effect upon the running. As a starting point we must then obtain 
the values of these couplings at a particular scale, for a chosen value of tan/3 and 
as{Mz)- Once we have selected these two parameters, we may determine the renor- 
malised masses mt^b,T{'<Tit) of the top, bottom and tau particles at the renormalisation 
scale mti by running the three loop QCD ® one loop QED RGEs through quark 
thresholds between and Msusy in the MS scheme [^. a5(/i < Mz) is actually 



run using a state of the art four loop QCD beta function ||3^. In fact, the four loop 
contribution only changes the predictions by a few parts per thousand. This fact has 
a limited significance because the coefficient functions required to extract as{Mz) 
from data are only known to at most three loops. The three gauge couplings are 
also run between Mz and in the MS scheme, assuming the particle spectrum of 
The Standard Model without the Higgs or top quark. The third-family MS Yukawa 



^mt denotes the running top mass in the DR renormaHsation scheme. 
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couplings may then be determined by p3 



\{.rnt) 



V sin P 

V cos P 



(37) 



where v = 246.22 GeV is the scale parameter of electroweak symmetry breaking. 
Above nit, we wish to use the RGEs for the MSSM contained in Appendix 1. However 
these are in the DR scheme and so at rrit we match all of the quantities obtained in 
the MS scheme to the DR scheme. A guess for the intermediate scale M/ is chosen 
and the gauge couplings and third family Yukawa couplings are run to this scale. 
Above Mi, the effect of the intermediate matter is felt and the RGEs change as 
prescribed in Appendix 1. The couplings are run up in scale until either it becomes 
non-perturbative (which we take to be greater than 4) or until (71 (yu) = (72 (/^)- In 
the first case the value of Mj chosen is abandoned and in the second a prediction 
for aslmt) is obtained by using the gauge unification condition. This is implemented 
by setting 5^3 (/x) = (72 (/^) and then running all of the couplings down to mt taking 
the intermediate matter into account and integrating it out of the effective theory at 
Mj. We may now iterate the above procedure using the previous predicted value of 
gsirrit) as an input each time until 03 (mt) converges and we have a value consistent 
with gauge unification with the intermediate matter at the guess value of Mj. The 
above procedure is then repeated for different values of Mj until a value is found in 
which the gauge couplings and unification scale satisfy Eq.|l|, i.e. the constraint of 
string scale gauge unification. It is a simple matter to re-convert 03 (m^) back into 
the MS scheme and run back down to determine as{Mz)- Thus, for a given tan/5, 
we now have the predictions Mj,as{Mz) that come from the assumption of string 
scale gauge unification. The conditions in Eqs.^,^ are then employed to remove 
any models which do not agree with the data or fit into the type of models of fermion 
masses being considered. 

Tables ^,|| display the results of the algorithm described above for tan /3 = 43, 5 
respectively. The constraints we impose upon the models are so tight that out of the 
tens of thousand models examined, only a few models pass the constraints in each 
case. Note that for tan/? = 43, only two of these give as{Mz) predictions within 
la of the central value. The minimum number of extra vector multiplets added is 
32. Varying mf sin^ 6'^ between their la limits can make a difference to as{Mz) 
predictions of ~ ±0.002. The Mj/Mx prediction is hardly affected by the change in 
the input parameters. Varying ai(M^) within its la limits makes only a negligible 
change to the predictions. We note that the greatest uncertainty in our two-loop 
calculation is likely to be that due to threshold effects. Until now we have assumed 
that all of the intermediate matter lies at one scale Mj. While this is in some sense 
the simplest scheme, in general there could be some splittings between the different 
types of additional matter. One may naively expect these to not span more than one 
order of magnitude, but even given this constraint there could be significant errors 
due to the non-degeneracy. Non-degeneracy of the superpartner spectrum could also 
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0.5579 
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0.1232 


0.5946 


0.5644 


5 


10 


10 


10 
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0.1231 


0.6098 


0.5716 



Table 3: Predictions of models that successfully unify the gauge couplings at Mx and 
provide enough intermediate matter to build a model of fermion masses for tan j3 = 43. 
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0.3737 



Table 4: Predictions of models that successfully unify the gauge couplings at Mx and 
provide enough intermediate matter to build a model of fermion masses for tan (3 = 5. 

cause errors in the predictions. If we allow the presence of sexton fields, there are 
some additional possible models that the algorithm just described will not find. These 
are models with equal one- loop beta functions above Mj. In this case, the one- loop 
algorithm fails because the solution to gi{fi) = (72 (yu) = fi'slyu) is not unique. Above 
Ml, the gauge couplings have slopes that differ by small two-loop effects. It is a 
simple matter to demonstrate that the models 

HQ = 0, ns = 24, n2-nD = 20; 

uq = I, ns = 16, n2-nD = H; 

uq = 2, ns = 8, n2-nD = 2; 

Uq = 3, ns = 0, n2-nD = 7; (38) 

have the property of equal one- loop beta functions above Mj. Some of these models 
were investigated with an accurate version of the two-loop algorithm. Table ^ displays 
the predictions of a subset of the successful models in Eq.pB|. 

4 Infra-Red Fixed Points 

We now turn to the question of infra-red fixed points (IRFPs) of the dimensionless 
Yukawa couplings for the class of models which are consistent with the generation of 
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Table 5: Predictions of models that successfully unify the gauge couplings at Mx and 
provide enough intermediate matter to build a model of fermion masses for tan f3 = 43. 
The models shown here belong to the special class of models that possess equal one- 
loop beta functions. 



acceptable textures via spaghetti mixing, and string gauge unification. Our discussion 



follows that of the IRFPs for the Ross model of Higgs mixing [|18|. The basic idea 
behind this approach is that when the large Yukawa couplings between heavy fields 
and Higgs' are renormalised, their low energy values may not be sensitive to the high 
energy ones, and so a rough prediction of the value may be made. This is the situation 
around an infra-red stable fixed point, where combinations of dimensionless couplings 
in the theory stop changing with renormalisation. Such fixed points are very welcome 
in this approach since our textures result from a large number of unknown Yukawa 
couplings, which would otherwise render this approach quite unpredictive. 

There is another reason for looking for infra-red fixed points in the models dis- 
cussed here. The presence of a gauged family symmetry such as [/(l)j5f is in principle 
quite dangerous since its presence can lead to large off-diagonal squark and slepton 
masses which can mediate flavour- changing processes at low energy. In particular the 
D term associated with U{l)x is in general only approximately fiat due to lifting 
by soft supersymmetry breaking terms, and this can lead to family-dependent squark 
and slepton masses with unacceptably large mass splittings. This is a generic problem 
of any model with a gauged family symmetry, however the U{l)x symmetry here is 
non-asymptotically free with a large beta function so that its gauge coupling rapidly 
becomes very small below the string scale, leading to small X gaugino masses. 



It has been suggested [|T^] that the possible infra-red structure of the theory could 
help by relating the soft scalar masses to the small gaugino masses, thereby making 
them naturally smaller than the squark and slepton masses, or by enforcing < 6 >=< 
^ > as an infra-red fixed point of the theory. We refer the reader to ref. |]l6l for more 
details. Here we shall only focus on the IRFPs of the dimensionless Yukawa couplings 
however. 



4.1 The Top Quark Yukawa Coupling 

The first step in finding the IRFPs of the theory is to construct the RGEs of the 
dimensionless Yukawa couplings of the theory. In supersymmetric theories this task 
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Figure 6: Contributions to the wavef unction renormalisation of the top quark Yukawa 
couphng. 

is made quite simple, at least at the one loop level, by the observation that only the 
wavefunction diagrams contribute to the RGEs. The vertex contributions vanish due 
to the non-renormalisation theorem. This allows the RGEs to be constructed in a 
very straightforward manner. To take a simple example, consider a toy theory which 
only involves the top quark Yukawa coupling in the superpotential: 

W = hQfH2. (39) 

Defining the Yukawa and three gauge coupling parameters as 



and the scale variable as 

t = -H^?) (41) 

we can write the RGE for the top quark Yukawa coupling as 



dt 



Y\Nq + Nt. + Nh,) (42) 



where Ni are the wavefunction renormalisation contributions from each of the three 
legs of the vertex. In the toy model the wavefunction diagrams are shown in Fig.^. 
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The wavefunction renormalisation contributions are explicitly, 

Nq = Y.'^C2i{Q)di-Y'' 

i 

i 

Nh, = Yl^C2i{H2)di-3Y'' (43) 

i 

where 2C2i{R) is the quadratic Casimir of the representation R under the i — th 
gauge group factor of the MSSM, arising from the gauge boson exchange corrections, 
and the multiphcity factors in front of the terms are due to doublet and colour 
counting for the particles going round the loop. Thus the RGE is explicitly 

rlY'^ 

—- = Y\Y1 ^C2i{R)d, - 6Y^) (44) 

Now let us assume that in our toy model all three gauge couplings were equal, and all 
three gauge beta functions were equal (quite unrealistic for the low energy Yukawa 
couphng, but typical of the situation near the string scale). Then the RGE may be 
written as 

HY^ 

^ = Y'^ira - ey'^) (45) 
at 

where we have defined 

r = J2'^C2i{R) (46) 

i,R 

where R runs over all fields involved in the coupling Y'^. We have written the three 
equal gauge couplings as d. We now write the one loop gauge running as 

'4-^ (4r) 

and define the ratio of Yukawa to gauge coupling as 

R^=— (48) 
d 

Then the RGE for this ratio is: 

dR^ 



^dR^{r + h- GR") (49) 



dt 



and the Pendleton- Ross fixed point is given by = 0. This condition can be 
achieved by 

r + b-6R'' = 0. (50) 
In terms of wavefunction renormalisation parameters, the RGE can be expressed as: 

d 

= R\Nq + Nt. + Nh, + dh). (51) 
22 



The fixed point condition can be expressed as^'': 

Nq + Nt. + Nh^ +ab = 0. (52) 

4.2 The General Model 

Ross [|r^ applied the above techniques to find the IRFPs of the Higgs mixing modeL 
We wish to extend the discussion to the more general situation of ng, n^/, n^, ri^, 
copies of the vector representations 

(g(^)+Q(-^)), (u^^y^+u'^-y^), (D^^y^+D'^-y^), (l^^^+l^-^)), {e^^^ + E'^-y^) (53) 

plus uh^ , copies of the vector representations 

{Hi^^ + Ht'^), {h['^ + Hi-'-"^) (54) 

where we have labeled the vector fields by their X charges x, z. For the special case 
of z = the Higgs are identified with the two MSSM Higgs doublets, 

Hf\ Hi'^ (55) 

and do not have vector conjugates. The above fields are in addition to the three chiral 
families of quarks and leptons which we label as: 

Q„ D'j, L„ £;| (56) 

where we label these fields by the family index subscript i,j = 1, ■ ■ - 3, but do not 
label their X charges (for i,j = 1,2,3 the X charges are -4,1,0, respectively, as 
discussed earlier). We introduce the X charge breaking singlet Higgs fields {9q + 
Oq), {9u + eu), {Od + Od), {Ol + ^l), {Oe + Oe) plus {O^ + h). (0, + ^i) which change 
the X charge of the particular field by 1 or -1. We also introduce MSSM singlet 
Higgs with X = 0: $q, ^u, ^l, plus $2, $1 whose VEVs are responsible for 
the heavy vector masses at a common scale Mj. 

The most general model is then defined by the gauge group, 

SU{3)c X SU{2)l X U{1)y x U{1)x (57) 

with the superpotential involving the chiral quarks and leptons containing the trilinear 

terms 

i,j id i,j 

^^It is worth noting that an alternative fixed point has recently been proposed by Jack and Jones 
[Q, in which = is achieved by the more stringent condit ions liVg = Nf: = = —\ah. 
These conditions may be expressed in a more general way, and are valid to two loops. However they 
are not satisfied for our simple toy model, and they will not be of use for our more general model. 
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i.y 



I'j 

+ E ^q.^qQ.Q^'"-'^ + E sc/,^c7f/;f/^(-^-^) + E SDfiDDp<-y-^^ 

i j 

+ E ^q.^qQ^Q^-^^'^ + E ^c/,^c/f/;f/'=^-^+'^ + E SD,eDDp'=^-y+'^ 

i j j 

+ E SLAL^L^~"'''^ + E SE,eEE]E<-y+'\ (58) 
* i 

In Eqs. |58| - |59| , it is to be understood that the X charges of the fields in each couphng 
must add to zero and that this decides the superscripts that are not summed over. 
This is true for all superpotentials and wave-function renormalisations listed in this 
paper. We neglect some terms in the superpotential that are not banned by the sym- 
metries we have listed so far. Some of these are undesirable in terms of reproducing 
the correct phenomenology, and so we appeal to the extra U(l) symmetries that tend 
to come with string-derived models to ban these terms. The remaining terms in the 
superpotential involving the extra vector states and Higgs are: 



W2 = Y.h-y)Q^''^u<y^H^^ + Y.h-y)Q^^^D<y^Ht 



(^) 

x,y x,y 

+ E h.y)L^''^E^^y^H[^^ + E h.)(y)Q^-''^U'^^-'^H, 



i-y) 

x,y x,y 



+ E k^)iy)Q^-''^D^^''^Ht'^ + E kx)iy)L^-^E^^"''>Ht^ 

x,y x,y 
z z 

+ E^QW-^Q^^'^O^""^ + Y.rjj(y)^uU<ym<-y^ + E^d(.)$d^'^'^5'^"'^ 

X y y 

+ Y.^Li.)^LL^^^n-^^ + Y.rE(y,^EE<y^E<-y^ 

X y 

z ^ z ^ 

2 ^ 2 ^ 

+ E ^qw^qQ^'^Q^"'"'^ + E ^q(^)^"qQ^^)q^-"+'^ 

X X 

+ E suiy^euu^^^y^u'^'^-y''') + E %(.)^i/f/'=^^)f/'=^^^+') 

y y 

+ E sniy)eDD<yW<-y-'^ + E soiy)enD<yW<-y+'^ 

y y 

X X 
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+ J2 SEiy)OEE<y^E'^^-y-'^ + J2 SEiy^eEE'^^y^E^-y^'^ (59) 

y y 

where in the first two hnes of Eq.|59|, X symmetry requires that z = —{x+y). Since the 
fields above are being labeled by their X charges, the limits of each of the summations 
will depend on the particular model under consideration. The family indices range 
from i,j = l,---3. However in specific models only a subset of the fields will be 
present, and consequently not all of the terms will be present. For the moment we 
prefer to keep the values of uq, nu, nD,nL, ue and nn^ , general, however. Also, we 
have not written the most general superpotential allowed under the gauge symmetry 
since 9e could couple to the vector quarks, for example. It is possible that 6'q,(7,d,l,£;,i,2 
are identified with just one superfield and that 4'q,u,d,l,e,h^_,H2 ^'^^ ^-Iso identified with 
one superfield (and similarly for the conjugate singlets). 

The one-loop RGEs for the couplings are 
dF& 

= 4(iVQ, + iV^c + iV^,.,+a6) 

dl^ 

= R\^{NL^+NE^^+N^,^.,+&h) (60) 
where we have assumed the gauge couplings are approximately equal so that 

a = -CKi = 012 = ct3. (61) 



The full list of wavefunction renormalisations are given in appendix 2. The fixed point 

-)h r)k T)l 



conditions for the chiral quark and lepton couplings, R^,,R^,,R\,, are listed below: 



NQ^ + Nuc + N^(.) + ah = 



'2 



Nq^ + Nd^, + N^i.) +ah = Q 



NL,+NE'^ + N^(.)+dih = (62) 



Rfx)jy R\x)jy th^t involve a mixture of chiral and vector fields. Also we require a 
fixed point for the couplings Rf^y), Rfxv)^ R\xv)y involving purely vector fields. Similar 



We also require a similar fixed point for the couplings -Rf(y), Ri(y)y R[{y)i ^'^^ R(x)j^ 

^{xy)-> ^^{xy)-> ^^{xy)-> 

fixed point conditions apply to the conjugate vector couplings, as well as all the sin- 
glet couplings. So every trilinear coupling will have a fixed point condition which is 
expressed in terms of the wavefunction renormalisations, similar to the above condi- 
tions. A fixed point is achieved when all the conditions are simultaneously satisfied. 
Note this assumes that none of the R^-^'^ couplings are zero at the fixed point, another 
set of possibilities allowed by Eq.|6y. We will not consider this here since many of 
the preceding arguments relied on the dimensionless couplings being ~ 0(1), rather 
than approximately zero. We merely note that in general there are 2" fixed points 
in this multi-dimensional system of n couplings, all but one of which involve some of 
the dimensionless couplings being zero. 
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4.3 Conditions for Infra-Red Attractiveness 



We now write the RG equations as 



dRijt) 
dt 



aRi{t) 



b)-j:S,,R,{t) 



(63) 



where Ri now denotes the ratio of any Yukawa couphng i to the gauge couphng 
(squared), as prescribed by Eq.^ We have written = 2J2C2{Rx), where the sum 
runs over simple gauge groups and the representations Rx under those gauge groups, 
X corresponding to the field that labels of Ri in Eq.^. The fixed point condition 
is then satisfied when the right hand side of Eq.^ is zero for all i. First, we assume 
that none of the Ri is equal to zero at the fixed point, in which case 



J2S^JR■=ri + b, 



(64) 



where we have denoted the value of Rj at the fixed point as R*. The problem of 
locating the fixed points becomes a straightforward problem in linear algebra, albeit 
involving a large number of dimensions, corresponding to the large number of trilinear 
Yukawa couplings. The fixed points are given in principle by inverting the matrix Sij, 



R* = j:iS'%irj+b) 



(65) 



To determine the infra-red stability of the system in Eq.|6^, we need to Taylor expand 
it around the fixed point given in Eq.^ We can then drop all except the linear terms, 
the resulting system of which allows an algebraic solution and can thus be tested for 
infra-red stability. We therefore make a change of variables to pi{t) = Ri{t) — R*. 
The RGE Eq.Bl then becomes 



dpijt) 
dt 



-a{t){p,{t) + R*) 



(66) 



where we have substituted the fixed point values of R* from Eq.^. When we drop the 
quadratic term in Eq.^ and change the independent variable from t to a by Eq.| 
we obtain the linearised system 

dpi{t) 



(67) 



Eq.^ then describes the behaviour of the trajectories as they approach the fixed 
point. It has solutions 

p,■(^)=x,(«(^))"^ (68) 

where Xj, \j satisfy the eigenvalue equation 



A.. = _R*Q. . 



(69) 
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Because the expanded RGE Eq.^ is linear, the general solution is a linear combina- 
tion of each pi in Eq.|67|. 

For 6 > as in these models, a decreases with decreasing renormalisation scale 
/i. For the fixed point to be infra-red stable, we require every eigenvalue Aj to have 
a positive real part, since then pj — >■ as decreases. This simply translates to the 
condition that every eigenvalue A of A must possess a positive real part. Complex 
eigenvalues always come in complex conjugate pairs, as do their associated eigenvec- 
tors. Writing Aj = kj + zsj, where fcj, Sj are real, the solution in this case is 



(70) 



Eq.|70| describes a spiral-like trajectory, the distance to the fixed point being controlled 
by a'^K Thus kj must be positive for the trajectory to be infra-red stable. 

Aj = corresponds to a direction in coupling space which is neither attracted 
nor repelled by the fixed point. For each of these directions there should be one free 
parameter in the solution to the fixed point equations. Thus, the dimension of the 
null-space of A gives the number of free parameters. The free parameters embody the 
information on where a solution lies along this direction (and are set by the initial 
boundary conditions). In the class of models presented here, this corresponds to some 
information about the string scale being retained at lower energies. In fact, it can 
be shown that the null-space of S is the null-space of A and so the number of zero 
eigenvalues of S fixes the number of free parameters. If the conditions for infra-red 
stability are not met, the fixed point is either a saddle point or an ultra-violet fixed 
point and so the fixed point will never be achieved at low energies. We will see in the 
following specific models, examples of infra-red stable and saddle point behaviour. 
We will also see that the null-space directions occur because of degeneracies in the 
fixed point equations. 



4.4 Example 1: Higgs Mixing Model 

As a first example of the general results, we calculate the fixed point solutions and 
the infra-red stability in a Higgs mixing model similar to that proposed by Ross [|18|. 
In this model there are nn^ = 10, = 10 copies of the vector representations 

(H^'^ + at^), {Hi'^ + Hi-^), (71) 

in the model, which means 77,2 = 20. Ross also included some colour triplets which 
served the purpose of increasing the gauge unification scale although not enough to 



be consistent alone with string-scale gauge unification. We saw in Eq.^ that such 
a model has string scale gauge unification if ns = 24, but we shall ignore the exotic 
sexton representations in the following analysis. The superpotential of the model is 
then 



i,j id 
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2 



2 = -2 

7 



(z)fj(-z-l) 



z=-2 
8 



(Z)fj{-Z + 1) 



z=-2 
7 



2 = -l 



z=-2 z=-l 

It is to be understood in the first three terms Eq.^ that 



2; = — Xcharge(z family) — Xcharge(j family). 



(72) 



(73) 



Therefore the Higgs which occur in the interactions with couplings [h, k, have 
charges as given below: 



( {h,k,l)uH['l ih,k,l)uHi'i (/i,A;,/)i3i/i5 \ 
{h,kJhM$ {h,k,l)22Hif^ {h,k,l)23H[-^''> 

1^ (/l,fc,/)3li/g {h,kJ),2Hi^'^ {h,kJUH^,2 J 



(74) 



The above Higgs having direct couplings are only a subset of the full list of required 
Higgses: 



o-(8) fr{-S) rrO) fr(-7) o-(6) o-(-6) ^.(5) o-(-5) o-(4) f.(-4) 
-"1,25-"!, 2 5-"l,25-"l,2 5-"l,25-"l,2 !-"l,2!-"l,2 5-"l,2)-"l,2 5 

o-(3) o-(-3) o-(2) o-(-2) „-(!) fr(-l) 7t(-1) ttW TT{--i) frCi) 
-"1,25-^1, 2 5 -"1,25 -"1,2 5 -"1,25 -"1,2 5 -"1,2 ) -"1,2? -"1,2 )-"l,2 



(75) 



The one-loop wavefunction renormalisations from Appendix 2 are: 

3 . 1 . 



Nr 



N, 



3 2 18 '-^ ^ * 



-0:3 H — 
3 ^ 9 



«i-2E^.? 



iVD^= = -d3 + -di-2^r, 
^L. = ^«2 + ^«i-E4 

Ne^ = 2di-2^FA 



(76) 



The Higgs wavefunction contributions are: 



"2 



l^2 + ld^-3YY,';-YY!^ 

3 . 1 . 
-ao H — «i 
2 2 
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(z) 



Ngi-.) = -a2 + -ai-Y "i -Y -Y 

: -02 + 2*^1 -Y "2 -Y "2 -Y "2 



The wavefunction contributions for the singlets are: 



z z 

-2j2Y'"f, Ne,=-2J2Y'" 



0Hl 



^ No,. 



(77) 



(7J 



Following Ross, we assume a general symmetric form for the matrices I, k, h at 
the string scale and use the fact that the RGEs Eq.^ respect this form. Like Ross, 
we also ignore all of the couplings in Eq.^ that do not involve MSSM states. The 



solutions of the fixed point equations Eq.^ applied to this model are 




-y - z 



103 _ 

72 c 
295 _|_ 

192 ~'~ ■ 64 
, 2^5 , ljU _ 
\ 192 ^ 64 "^"^ 



295 I ll_b _ o 

192 ~'~ 64 y 

295 I 11_6 _ o 
192 64 "-^ 



1297 _|_ 
864 



X 

z 

5 b 
32 



■ Z 



295 _|_ 11 b 

192 ~'~ 64 
295 _|_ 11 b 

192 ~'~ 64 



103 
" 72 




Eq.l79| shows that there are 3 undetermined parameters at the fixed point. Note that 



if any of the entries of the matrices in Eq.|79| are negative, the coupling corresponding 
to the negative solution will tend to zero. This could violate the assumptions about 
the fundamental dimensionless couplings all being of order unity and hence the choice 
of X charges required to reproduce the phenomenology. We therefore demand that 
x,y,z must satisfy the constraints 



{x + y,x + z,y + x} > 
0<{x,y,z} < 



103 b 

216 24 
295 116 

576 ^ 192' 



^0) 



When we check the solution for infra-red stability we find that the eigenvalues of A are 
positive for solutions that satisfy the above conditions. Thus, the fixed point solution 
identified is completely infra-red stable, with three undetermined free parameters. 

So far, in the limit that singlets are ignored, our results are in agreement with those 
of ref. |T^. Although the stability question was not explicitly addressed [|I^ we find 
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that the fixed point is stable in the infra-red hmit so all is well. Now we must consider 
the effect of the singlets. In the Ross model the result was quoted that the singlet 



couplings of the Higgs Mixing Model are approximately flavour independent. We flnd 
that this is only valid if the extra Higgs doublets which do not have direct couplings 
to fermions are ignored. To be explicit, Ross considered a model with the only extra 
Higgs states being 

tt{S) fri-S) ttW fr(-4) rr(3) fr(-3) rr{2) rr(-2) ttW fr(-l) /oi \ 

-"l,2;-f^l,2 ) -"1,2; -f^l,2 ) 1,2 ; -^^ 1,2 ) 1,2 ; -^^ 1,2 ) 1,2 ; -^^ 1,2 • \.°^J 

However we saw earlier that the full list of Higgs states in Eq.^ is required for correct 
Cabbibo mixing and CKM mixing. The full Higgs mixing model is analysed in Ap- 
pendix 3 where we solve the 80 simultaneous equations for 80 unknowns (keeping the 
matrices k, I, h symmetric) to determine the predictive properties of the model. The 
solution detailed in Appendix 3 shows that the solution has 27 undetermined param- 
eters. We now discard the model because we have more unconstrained parameters 
than data points on fermion masses and mixings. 



4.5 Example 2: Quark-Line Mixing Model 

We next consider a model with mixing along the Q doublet line provided by the extra 
flelds 

Q'^'\Q^'\Q^'\ g(^), g(°), Q^-'\ Q^-'\ Q^''\Q^~'^ (82) 

and their partners in the vector-like pair. Such a model by itself is not expected to be 
realistic since it does not account for lepton masses, but it may be regarded as part 
of a fuller model such as the nq = 9,nu = 4:,nD = 10, ue = S example in Tables ||,^ 

The superpotential of this example is explicitly: 

x=—A 

+ E SQ,.,eQQ^^^Q^-^-'^ + E -SQi^^eQQ^-'^Q^-^^'^ (83) 

x=— 4 x=— 3 

The wavefunction renormalisations for the chiral quarks and leptons and MSSM 
Higgs doublets are (see Appendix 2): 

Nq, = ^cia + ^«2 + - ^3^3 - r^«3 - F^-Q3 

8 3 1 
Qi 3 2 18 
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— 




Nuf 


= 




Nh, 


— 




Nqw 


= 




Nq(-i) 


— 






— 




NqW 


— 


Nq(. 


:ytl,0,-l,4) 


— 






= 




Nq(-2) 


— 




Nq(-i) 


= 




Nq(o) 


— 










-2,-1,0,1) 























7T«3 + 7T«i - 2r3^3 - 2Y(>X 



33 ^-^ (0)3 

8 8 

3^3 + g«i - 2y('^i)2 

8 8 

^d3 + + - F(^)3 - ^'=(0) - r^'^(o) - F^"'5(o) 

^«3 + + - Y('L,), - y'^(-) - y^'^(-) - 
Ids + + ^d, - yH^, - - F^"«(^) 

^«3 + + ^di - Y'^'m - y^'^d) 

1^3 + + ^di - y^'^w - y^'^w - y^'^w 

3 ^ 2 18 

-^3 + -d2 + —di - Y'q^'^ - y'«(^) - Y'q^'^ - Y'^2 
3 2 18 

-a~3 + -^2 + — di - yod) - y^"Q(2) - y^Q(°) - y-Q3 

3 2 18 

-da + -d2 + — di - yo^"^ - y"o(^) - y«'-^) - 

3 2 18 

-da + -d2 + — di - y'q^-'^ - y'Q(°> - y^o'"^' - y^'^a 

3 2 18 

^«3 + ^d2 + ^di - yo(--) - y^"Q(-+i) - y^Q(— 1) 

-6 j2 

a;=-4 
3 

-6 ^ y'«'-) -6(y«i +y"«2 +y^«3), 

x=-4 
4 

-6 ^ y'«^") -6(y«2 +yo3). (84) 

a;=— 3 

The fixed point conditions for the couphngs Rz^fi%, 
, are: 

NQ,+Nuc + NH^+dh = (85) 

Nqw + iVt/c + ATj^^ + = (86) 

NQ,-r, + Nuc + Nh, + db = (87) 

Nqw + Nui + Nh, + db = (88) 
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^1(3) + ah 



^e-Q + Nq^ + A^Qfo) + ah 
% + A^Qg + Nqw + a& 

+ ah 
+ (56 








0, 
0, 



X e {-4, 

X G {3, . . 
0, X G {4, . . 



..,4} 

,-4} 
,-3} 



(89) 
(90) 
(91) 
(92) 
(93) 
(94) 
(95) 
(96) 



If we were to ignore the contribution of the singlet sector then the fixed point 
equations for the couphngs R^^jR'^j, Eqs.^ and lead to the matrix equation: 



/ 6 5 3 3 \ f 

5 6 3 3 

3 3 6 3 

V 3 3 3 6 y 



-^33 

^(0)3 

^(-1)2 
vh 
1 1 



\ 



V ^(4)1 J 



1 
1 

V 1 J 



r 



QUH2 



+ h)a 



(97) 



where r'^^^^ = 88/9, and we have assumed all the gauge couplings to be equal. Upon 
inverting the matrix we find, 



vh 
-'^33 

vh 

^(0)3 

vh 

^(-1)2 
vh 
^(4)1 



3 
5 

V 5 / 



.QUH2 



+ 6)_ 



63 



-a 



(9^ 



for the fixed point solutions of the couplings. Note that this solution has a global 
SU (2) flavour symmetry in the Yukawa couplings of the two lightest families, unlike 
the Higgs mixing model for example [|l^ . The reason that it is present in this model is 
that there is a single Higgs doublet which is common to all the fixed point equations, 
as compared to the Higgs mixing model where a different Higgs couples in each entry 
of the Yukawa matrix. When the singlets are included they will explicitly break 
the global SU{2) flavour symmetry, as we discuss below. Note that in some recent 
models, such a symmetry is assumed as a starting point We then checked 

that the system of RGEs in Eq.|62| is infra-red stable in this case by determining the 
eigenvalues of the matrix Aij. These come out to be 6 + r'^^^^, {h + r'^^^2)/21, 
5{b + r'^^^^)/21, 5{b + r'^^^^)/21, so the fixed point is encouragingly infra-red stable 
in all four independent directions. 



Once the singlets are included the above fixed point in Eq.pB| are be modified. If we 
return to Eqs. ^5| - |96| we see that there are the same number of equations as unknowns, 
so the whole system may be regarded as a large matrix which may be inverted along 
the above lines. From Eqs. |85| - p6| the following relations may be obtained. 



Nr 



N 



Nq, 



N 



Nr 



N 



(99) 
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= i^i^Ne^+N-e^) (100) 

= Nq, - Nq, = Nu^ - Nui 

= ^(Nq, - Nq,) = ^{Nug - Nu^) 

= liNQ, - Nq,) = ^{Nu^ - Nu^^) (101) 

These relations are formally quite model-independent: they apply to any model with 
hq = 9 provided Eq.^ holds, regardless of the number of additional states. However 
the implications of these relations will depend on the particular model under consid- 
eration since the wavefunction renormalisations have model dependence. For instance 
in this particular example, we can immediately see that the previously obtained fixed 
point based on ignoring the effect of the singlets is not consistent with these equalities. 
For example it would imply {Nu--Nu-) oc (5-5) = and {Nu--Nu-) oc (9-10) 7^ 0, 
although the two relations are approximately consistent. 

In Appendix 4 we give the fixed point of this example, including the singlets. 
The fixed point solution has 9 undetermined parameters and so again, we see the 
predictivity of the model is severely lowered by the inclusion of the singlet couplings. 
This behaviour is rather similar to that encountered in the Higgs mixing model, and 
we therefore expect that it may be a general feature of models of this kind, once 
the singlets are included. There is also another serious problem with the fixed point 
solution, given that it is impossible to pick values of the 9 free parameters that predict 
all of the constrained couplings to be positive. Thus, the fixed point identified cannot 
be realised by nature as one or more couplings will be forced toward zero in any 
infra-red stable solution. 



5 Conclusions 

We have explored a scenario in which the minimal supersymmetric standard model 
(MSSM) is valid up to an energy scale of ~ 10^^ GeV, but that above this scale the 
theory is supplemented by extra vector-like representations of the gauge group, plus 
a gauged U{l)x family symmetry. The basic idea of our approach is that the extra 
heavy matter above the scale ~ 10^^ GeV may be used in two different ways: (1) to 
allow (two-loop) gauge coupling unification at the string scale; (2) to mix with quarks, 
leptons and Higgs fields via spaghetti diagrams and so lead to phenomenologically 
acceptable Yukawa textures. 

We considered models in which there are enough heavy vector representations 
to give every effective MSSM-type Yukawa coupling a non-zero value. Using this 
constraint (detailed by Eqs.^^^), plus the further condition that the mass scale must 
not be too far below the string scale, we performed a two-loop string gauge unification 
analysis which yielded 8 models that satisfy these conditions for tan /? = 43 and 2 for 
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tan /3 = 5 in Tables ^ and ^. For example ng = 9,nu = 4:,ni:) = 10, ue = 8 (all other 
Hi zero) satisfies the constraints of string unification independently of tan /?, and also 
has enough heavy matter to enable Yukawa textures to be generated via spaghetti 
diagrams. An example of a different solution is n2 = 20, ns = 24 where n2 = 20 may 
be interpreted as being due to uhi = 10, = 10 as required in the Higgs mixing 
model. 

Because the dimensionless couplings are of order 1 and because the RG running 
of the gauge couplings above Mj is steep, one might hope the dimensionless couplings 
would be forced toward numerical values approximating a fixed point at Mj. This 
would allow us to make numerical estimates of the values of the fermion masses and 
mixings at low energy. Motivated by these considerations we constructed the super- 
potential for a general model involving intermediate matter, and various Standard 
Model singlets that provide the U(l)x breaking. We obtained the one-loop RGBs 
for the general model, and then obtained conditions for stability of the fixed point, 
showing that the direction of stability in terms of the renormalisation scale depended 
on the eigenvalues of the coefficient matrix of the fixed point equations. 

Having discussed the general case, we then investigated the fixed point of two 
examples in some detail: a Higgs-line mixing model with = 10, = 10 and 
a quark-line mixing model with nq = 9. Both models have infra-red stable fixed 
points in the approximation that the couplings involving the singlets are ignored. 
However when the singlets are included in the analysis we found that the number of 
undetermined parameters grows (from 3 to 27 in the Higgs-line mixing model, and 
from to 9 in the quark-line mixing model). The spirit of the fixed point approach 
is that the model's couplings will focus down on a particular set of values at a lower 
energy scale independent of the starting conditions at a higher energy scale. In the 
full Higgs-line mixing model, the predictive fixed point solution is only valid given 
certain constraints on the values of the couplings at high energy, and so this seems 
to be in opposition to the spirit of the fixed point approach. In the full quark-line 
mixing model, the fixed point solution demands a negative value for at least one of 
the dimensionless coupling (squared) Ri. Any Ri that have negative fixed point so- 
lutions will actually tend toward zero, violating the initial assumption of order one 
dimensionless couplings that led to the construction of the U{l)x model. If the full 
Higgs-mixing model is to have a solution for which the couplings squared are all pos- 
itive, a complicated set of constraints must hold upon the boundary conditions at 
the higher energy scale. While the examples investigated here predict rough orders 
of magnitude for the fermion masses and mixings, it appears that detailed attempts 
to make them predictive by the analysis of infra-red fixed points is problematic as 
the above cases illustrate. It is unclear how to avoid these pitfalls in the construc- 
tion of predictive models without the inclusion of some extra symmetry to increase 
predictivity, and/or boundary conditions provided by an explicit string model. 

The lesson learned is clear: it is not in general appropriate to ignore the singlet 
couplings which must be incorporated fully into the analysis. The most predictive 
scenario would be one in which the Yukawa couplings depend only upon a in the fixed 
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point solution. It is not yet clear how one could pick a model in which this is likely to 
be true, or how one could pick a model that possesses a completely infra-red stable 
fixed point with all Yukawa couplings non-zero. However it is possible that such a 
model is contained in the subset of the models in Table ^ which have not all been 
analysed in detail because of their individual algebraic complexity. The search for a 
completely realistic model, and the detailed comparison of its low energy predictions 
to data, is beyond the scope of this paper. 

The idea of being able to predict the entire fermion mass spectrum in terms of 
one or two free parameters is an exciting prospect, and we hope that the general 
results of the present paper will be helpful in this endeavour. While we were not 
able to make the models quantitatively predictive in the fermion mass and mixing 
sector, they simultaneously provide an explanation for string gauge unification and 
the hierarchies in the fermion mass sector. 
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Appendix 1: Two-Loop Renormalisation of MSSM Plus In- 
termediate States 



We derived the following RGEs for the two-loop evolution of the gauge couplings 
and third family Yukawa couplings in the case of additional intermediate matter from 
ref. [^. Note that we have neglected other Yukawa couplings as an approximation. 
This a good approximation for the bulk of the running which is between Mj and 
1 TeV, where the intermediate states have been integrated out of the effective field 
theory and the effective Yukawa couplings are all much less than 1 apart for \t,b,T- As 
we are not considering neutrino masses in this analysis, we assume that there are no 
neutrino Yukawa couplings in the effective field theory being considered. For now we 
must assume there are no extra couplings between the superfields of the MSSM and 
the extra matter for reasons of simplicity and computing time. Naively one might 
expect these couplings to only change the results slightly because they decouple after 
less than 2 orders of magnitude in renormalisation running and because they only 
affect the running of the gauge couplings at the two-loop level. Nonetheless, it should 
be borne in mind that these couplings could influence the results, particularly in view 
of the fact that these dimensionless couplings are expected to be of order 1. The 
equations are valid in the DR scheme. 
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where t' = In /i and /i is the DR renormahsation scale. 



We now detail the matching conditions between the DR and the MS schemes [p6|: 
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(103) 



where the superscripts denote what scheme the quantity is evaluated in. To a good 
approximation, g^^ = gf^^ and A:^^ = X^^'^ . Our running value of rrit is determined 
by 
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where m\^^^ is the value extracted from experiment. For m\ 
central values of as{Mz), we obtain rritimt) = 166 GeV. 
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Appendix 2: One- loop wavefunction renormalisation of the 
general model 



Here we give the wavefunction renormalisation contributions to the RGEs for the 
general model in Eq.^ In the following equations for the wavefunction renormali- 
sation of Ni, all sums are intended to be over couplings that multiply the field i in 
Eq.|5|: 
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i X 

The wavefunction renormalisations for the vector states are: 
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The wavefunction contributions for the singlets are: 
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Appendix 3: Fixed point solution of the Higgs Mixing model 
including singlets 

For the Higgs Mixing model corresponding to the superpotential in Eq.^ with 
symmetric inputs for the /, /i, k matrices at a high scale but including the singlet 
couplings, we now present the solutions to the fixed point equations: 
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Appendix 4: Fixed point solution of Quark-line mixing model 
including singlets 
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On solving the fixed point Eqs. |85| - p6| including the singlet sector, we obtain 
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(109) 



Eq. |109| shows that in fact 25 out of the 34 variables are constrained. This must 
mean that within Eqs. |85| - |96| , 9 of the 34 constraints exhibit degeneracy. Note that 
for R^q'^^^ > 0, we must pick i^^Q'"*' to be negative. 
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